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I. INTRODUCTION 



Recently, supersymmetry (SUSY) and shape invariance in quantum mechanical matrix mod- 
els have attracted much attention in the literature, e.g., Refs. [1-14] in various physical con- 
texts. To avoid confusion, we note that we here mean SUSY between two matrix Schrodinger 
operators intertwined by a matrix linear differential operator but not SUSY between two 
scalar Schrodinger operators in a matrix superHamiltonian intertwined by a scalar linear 
differential operator. 

On the other hand, the framework of jV-fold SUSY [15-17] has been shown to be quite 
fruitful among several generalizations of ordinary SUSY especially since the establishment 
of its equivalence with weak quasi-solvability in Ref. [16], for a review see Ref. [18]. Due 
to the facts that the Af = 1 case corresponds to ordinary SUSY and that shape invariance 
automatically implies weak quasi-solvability, A/"-fold SUSY contains both ordinary SUSY 
and shape invariance as its particular cases. 

Hence, it is quite natural to ask whether a formulation of A/"-fold SUSY is possible for 
matrix systems. To the best of our knowledge, there was so far only one such an attempt 
which corresponds to the Af = 2 case for 2 x 2 matrix models [4]. However, the analysis 
there was quite restrictive, was devoted mostly to the particular cases where the systems were 
available by two successive SUSY transformations, and resorted to the involved assumptions 
and ansatz. As a consequence, its formulation does not have a form which is suitable for 
discussing general aspects, especially those which were established and appreciated later 
after that work. 

In this article, we formulate for the first time A/"-fold SUSY for a system composed of 
matrix Schrodinger operators for all positive integral Af in such a general fashion that the 
recent crucial developments in the field are fully incorporated in the formalism. To see such 
a system actually exists, we construct as an illustration general 2x2 Hermitian matrix 2-fold 
SUSY systems without any assumption or ansatz. We find that there are two inequivalent 
systems, both of which are characterized by two arbitrary real scalar functions. Intriguingly, 
one of the systems does not admit reduction to the most general 2-fold SUSY scalar system. 

We organize this article as follows. In the next section, we generically define A/"-fold SUSY 
in quantum mechanical matrix models. Then, we investigate in detail 2x2 Hermitian matrix 
systems for Af = 2 in Section III. We explicitly solve all the conditions for 2-fold SUSY 
to obtain general form of the latter systems. In the last section, we refer to several future 
issues to be followed after this work. 



II. GENERAL SETTING 

A quantum mechanical system we shall consider here is a pair of n x n matrix Schrodinger 
operators 
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where I n is an n x n unit matrix and each V ± (q) is an n x n matrix- valued complex function. 
Let us introduce a pair ofnxn matrix linear differential operators of order Af 



d k 



k=0 



dq k 



AN ^-1 Ak 

p M=(-^in^+j:(-i) k ^- kW kWh 

k=0 



dq k 



(2a) 
(2b) 



where Wk(q) (k = 0, . . . ,Af — 1) are nx n matrix-valued complex functions. Then, the 
system (1) is said to be Af-fold super symmetric with respect to (2) if the following relations 
are all satisfied: 



P%H* - H±P% = 0, 



^N 1 N 



Pf rP ± = 2 X 



Af-l 



Af-k-l 



k=l 



(3a) 
(3b) 



where (k — 0, . . . , Af — 1) are n x n constant matrices. It is evident that in the case of 
n — 1 the above definition of A/"-fold SUSY reduces to the ordinary one for a pair of scalar 
Schrodinger operators. As in the scalar case, the first relation (3a) immediately implies 
almost isospectrality of H ± and weak quasi-solvability H ± ker Pj^ C ker P^. We note that 
in contrast to the formulation in Ref. [4] where its focus was on hidden symmetry operators 
characterized as the deviation from the second algebraic relation (3b) we have treated it as 
a part of the definition of A/"-fold SUSY. See also the discussion in Section IV, item 1 for its 
relevance. 



III. TWO-BY-TWO HERMITIAN 2-FOLD SUPERSYMMETRY 

As an example, let us construct the most general 2x2 Hermitian matrix 2-fold SUSY 
systems which are defined in the Hilbert space L 2 of two-component functions equipped with 
the inner product defined by 

(0,^)= [ dqtfiqMq), <f>,^ G L 2 (S), Set, (4) 
Js 

where the Hermitian conjugate f is as usual the combination of complex conjugate and 
transposition. The most general form of a pair of 2 x 2 Hermitian matrix Schrodinger 
operators is given by 

q /i=0 

where <Jq = I2 is a 2 x 2 unit matrix, <jj (i = 1,2,3) are the Pauli matrices, and V^(q) 
(/j, = 0, . . . , 3) are all real scalar functions. Components of 2 x 2 matrix 2-fold supercharges 
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have the following form 

d 2 ( 3 Ad 3 
p 2 = h~7^ + W *MW )tz + Yl W M°^ ( 6a ) 

d 2 d ( 3 \ 3 

P2+ = h ck? ~~ dq 1 ^ Wl ^ a » I + Yl w oM°i» ( 6b ) 

\/i=0 / ^=0 

where Wi^(q) and w ^(q) (/i = 0,...,3) are all real scalar functions. To investigate the 
TV-fold SUSY condition (3) for the M = 2 case under consideration, we first note that 
P 2 ~H~ — H +P 2 = implies P 2 + H + — H ~P 2 = and vice versa since they are Hermitian 
conjugate with each other with respect to the inner product (4), thanks to the choices (5) 
and (6). Hence, it is sufficient to study only the former. A direct calculation shows that the 
intertwining relation P 2 H~ — H + P 2 = holds if and only if the following set of conditions 
is satisfied: 

v; - v- = (?) 

3 

w'[ + 2w' 00 + 4V -' - 2 ^(V+ - V » ) = °> (8) 

< + 2w' 0i + AVr' - 2w u (V Q + - V ~) - 2w w {V l + - V~) = 0, (9) 

3 

E^«(V + ni = o. (10) 

j,k=l 

3 

+ 2v -" + 2 K V - w ^ V t ~ = °» ( n ) 

< + 2V~" + 2w u V -' + 2w w Vr' - 2w 0i (V+ - V ~) - 2w 00 (V+ - V~) = 0, (12) 

3 

E ^ kV + w o 3 (V k + + V k ~)] = 0. (13) 

j,k=l 

On the other hand, we find that the 2-fold superalgebra P 2 T P 2 = 4 [(H^ + C ) 2 + C x ] with 
Hermitian constant matrices Ct = Ylli=o Ck^c*^ (Cfc^ £ R, A: = 0, 1) holds for the upper sign 
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if and only if 

3 

W + = 3w' w - 2w 00 + J2( w ^) 2 - 4 ^oo, (14) 
AV+ = 3w' u - 2w 0i + 2w w w u - 4C 0t , (15) 

3 

E € ijk W lj( W lk ~ W 0k) = 0, (16) 
j,k=l 

3 

2V+" - + Co,) 2 - 4C 10 = 

3 

w'w - w oo + [ Wl » w % + w 'i, w ^ - wi,w' 0fl - (wo,,) 2 ] , (17) 

2V+" - 8(V+ + C 00 )(V; + + C 0i ) - AC U = < - < 

+ wio^ii + tyio^ii - wd wii + woow'u + w' w w 0i - w w w' 0i - 2w 00 w 0i , (18) 

3 

^ £ijk{w"jWik + w'^wok + w^w'm) = 0, (19) 
j,k=i 

and for the lower sign if and only if 

3 

W ~ = -w[ - 2w 00 + J2^ 2 ~ AC ^ ( 2 °) 
AVr = -w' u - 2w 0i + 2-wwWu - 4C 0i , (21) 

3 

^ tijkWijWok = 0, (22) 
j,k=i 

3 3 

2 V - 4 E(^" + C °») 2 - AC ™ = ~<o + E + ^< - (™o„) 2 ] , (23) 

2VT" - 8(l/ - + C m ){V~ + Co,) - AC U = 

- w'w + w' 00 wu + wq Q w' u + w' 10 w 0i + w w w' 0i - 2w ow i, (24) 

3 

E e ijk{w'ijWOk + WljKk) = °- ( 25 ) 
j,k=l 

The formulas (14), (15), (20), and (21) determine the form of the potentials and are 
compatible with the conditions (7), (8), and (9). Then, the conditions (10), (16), and (22) 
are identical with 

3 3 3 

E ^ijkWljWok = E e ijkWljW' lk = ^ ^jk^ljCok = 0. (26) 
j,k=l j,k=l j,k=l 

The most general solutions to the latter set of conditions are given by 

wu = C M Vi, w 0i = C 0i v , (27) 
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where v\ and vq are at present arbitrary scalar functions. Substituting (27) into (14), (15), 



(20), and (21), we obtain 

4Vq + = 3w[ - 2w 00 + (w w ) 2 + C\ Vl ) 2 - 4C 00 , (28a) 

4V+ = C 0l (Sv[ - 2v + 2w loVl - 4) , (28b) 

4V ~ = -w' w - 2w 00 + (w 10 ) 2 + C 2 {v x ) 2 - 4C 00 , (28c) 

AV~ = C 0i (-v[ - 2v + 2w loVl - 4), (28d) 



where C 2 = Eti(Coi) 2 - The solut ions (27) automatically satisfy (13), (19), and (25). With 
the substitution of (27) and (28) into the remaining conditions, (11) and (12) are respectively 
identical with 

w'w ~ w w w'( - 2(w' w ) 2 + Aw' 10 w 00 + 2w 10 w' 0Q - 2(w 10 ) 2 w' w 

- C 2 [vW{ + 2K) 2 - Hvo - 2^1^ + 2KoK) 2 + 4w l0 v lV [] = 0, (29) 
v 'i ~ w io v i - ^ w w v i - w w v i + 4w' 10 v + 2u) 10 v' + 2w' 00 v 1 + Aw m v[ 

- Aw 10 w' wVl - 2(w w ) 2 v[ - 2C 2 {v 1 ) 2 v[ = 0, (30) 

(17) and (18) are respectively with 

2w'" - 5(w' w ) 2 + 8w[ w 00 + Aw 10 w' 00 - 6(w w ) 2 w' w + A(w 10 ) 2 w 00 - (w 10 ) A 

- C 2 [5(v[) 2 - 8v[v - Aviv' + Qw[ (vi) 2 + 12w w viv[ - 8w 10 viV 

- Aw 00 ( Vl ) 2 + Q(w w ) 2 ( Vl ) 2 } - C 4 ^) 4 - 16C 10 = 0, (31) 
C 0i [v'l' - 5w' w v[ + Aw' w v + 2w w v' + 2w' 00 v 1 + Aw 00 v[ - 6wi K v i - 3{ w wY v 'i 

+ 2(w 10 ) 2 v + 4w 10 w 00 v 1 - 2(w w ) 3 v 1 - C 2 (^) 2 (3^ - 2v + 2w 10 v 1 )] - 8C U = 0, (32) 

and (23) and (24) are respectively with 

2w w - 4:W 10 w" - 3(w' w ) 2 + 8w' w w o + 4:W W w' 00 - 2(ww) 2 w' w - 4(wi ) 2 u>oo 
+ (w w ) 4 - C 2 [Av^'l + 3(v[) 2 - 8v[v - A Vl v' + 2w' 10 ( Vl ) 2 + Aw^v^ 
+ 8w w v lVo + 4w oK) 2 - 6(w 10 ) 2 K) 2 ] + C 4 K) 4 + 16C 10 = 0, (33) 

C 0i [v'l' - 2w'{ a v 1 - 3w' 10 v[ - 2w w v" + 4w' l0 v + 2wi v' + 2w' 00 vi + 4w ov[ 
- 2w w w[ v 1 - (w w ) 2 v' 1 - 2(w w ) 2 v - 4w w w 00 vi + 2(w w ) 3 v 1 

+ C 2 ( Vl ) 2 (-v[ - 2v + 2w wVl )] + 8C U = 0. (34) 

It is evident that the equations (32) and (34) have the trivial solutions Cu = Co, = 0. On 
the other hand, for non-trivial solutions the combination Cu/Coi :— C {i = 1,2,3) should 
not depend on the index i. 

Let us first consider the set of conditions (29) and (30). We find that the two combinations 
wi x (29)+C 2 t>iX (30) and v i x (29)+u>i x (30) are total differentials and thus are integrated 
respectively as 

2w w w'{ - (w' w ) 2 - 2(w 10 ) 2 w' w + A(w w ) 2 w 00 - (w w ) 4 

+ C 2 [2v lV '( - K) 2 - 2u^ (t; 1 ) 2 - Aw^v^ + 8w 10 vit;o 

+ Awooiv,) 2 - 6(^ 10 ) 2 (t;i) 2 ] - C\ Vl Y = D 1: (35) 
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and 

w" vi - w' 10 v[ + w 10 v" - 2w w w' w v 1 - (w w ) 2 v' 1 + 2(w 10 ) 2 v + 4w w w oVi 
- 2(w w fv 1 - C 2 (^) 2 (^ - 2v + 2w w v 1 ) = D 2 , (36) 

where D\ and D 2 are integral constants. It is easily checked that (35) and (36) are compatible 
with all the remaining conditions (31)-(34) if and only if 

D x = 16C 10 , D 2 = 8C U /C 0l = 8(7. (37) 

Hence, the only remaining problem is to analyze (35) and (36). They can be regarded as 
simultaneous linear equations for woo and t> . For the non-degenerate case vi ^ C~ 1 wio, 
they are uniquely solved as 

4 [(w m ) 2 - C 2 ( Vl ) 2 ] 2 w 00 = Ko) 2 [-2w 10 w'; + K ) 2 + 2{w w ) 2 w' w + {w w f + 16C 10 ] 
+ C 2 [2w w w'{o{vif + (w[o) 2 (vi) 2 - Awww'wVivi + 2(w 10 ) 2 v 1 v'I + (w 10 ) 2 (v[) 2 

- 4( Wl0 ) 2 ^ (t; 1 ) 2 - (w 10 ) V) 2 - 32^^! + 16C 10 (^i) 2 ] 

- C 4 K) 2 ptontf - K) 2 - 2^o(^i) 2 + («>io) 2 M 2 ] + ^ 6 K) 6 , (38) 

and 

2 [(wio) 2 - C 2 (wi) 2 ] 2 w = w 10 [w 10 w" vi - (wio) 2 ^! + Wiowi ^ - (w w ) 2 v" + (wio)M 
+ (iuio)V + 8(5^10 - 16Ci fi] + CV [-Wio(fi) 2 + ty'io^i^i + 2w io^i - ^ioK) 2 

- 2(w 10 ) 2 v 1 v' 1 - 2(w 10 ) 3 (v 1 ) 2 + SCv 1 ] + C 4 ( Vl ) 4 (v[ + w 10 v 1 ). (39) 

Finally, the general form of a 2 x 2 Hermitian matrix 2-fold SUSY system for the non- 
degenerate case is given by 

H + = - ~ + \ [3u/ 10 - 2w 00 + (w 10 ) 2 + C 2 ( Vl ) 2 } 

1 3 
+ -(3v[ - 2v + 2w 10 v 1 ) C 0i<Ji - C , (40) 

i=\ 

H ~=-\^2+\ - 2^oo + (w 10 ) 2 + C 2 K) 2 ] 

1 3 
+ -{-v[- 2vo + 2w 10 V!) C u°i - Co, (41) 

i=i 

d 2 / 3 \ d 3 

P2 = dg 2 " + y 10 + Vl ^ Coi<Ti J + W °° + V ° S C ° iCTi - ^ 

The two functions woo and t>o in the above can be eliminated by using (38) and (39). Hence, 
the system can be expressed solely in terms of the two functions w w and V\. 

For the degenerate case V\ = C~ 1 Wio, the two equations (35) and (36) are not independent 
and are equivalent with the following single equation 

4(u>io) Voo + Cvo) = -2w 10 w'( + K ) 2 + 4(w 10 ) V 10 + 4(u> 10 ) 4 + 8C 10 , (43) 
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with Cio = CC. Hence, we can again eliminate two of the four functions, e.g., v\ and vq- 
The general form of a 2 x 2 Hermitian matrix 2-fold SUSY system for the degenerate case 
is given by 



H + = 



^d7 + i^°" 2woo+2(Wio)2 ] 



i 

AC 



w 



Ko) 2 



w' w + 2w o + — 



ACC 



w 10 2(w 10 ) 2 (w 10 y 



i=l 



H~ = 



1 

AC 



3w' 10 + 2u>oo + 



w 



(Ko) 2 



10 



ACC 



iuio 2(wi ) 2 (wio) s 



^] Coi°"i — Co, 



A 2 

dg 



2 + Wio 



0i a i 



d^ 



+ w 00 + 



C 



! / \2 ^10 

w 10 - woo + {wio) ~ ~ h 



10 



/ \2 
10 ) 



+ 



2CC 



2w 10 4(wio) 2 {wio) 



i=i 



(44) 



(45) 



(46) 



It is interesting to note that in the limit — >■ (i — 1, 2, 3), the non-degenerate system 
(40)-(42) reduces to the most general 2-fold SUSY scalar system [16, 19, 20] while the 
degenerate system (44)-(46) does not. 



IV. DISCUSSION AND SUMMARY 

In this article, we have for the first time formulated generically A/"-fold SUSY in quantum 
mechanical matrix models and constructed the general 2x2 Hermitian matrix 2-fold SUSY 
systems without recourse to any assumption or ansatz. In addition to the detailed studies 
for larger n x n matrices (n > 2) and Af > 2 cases, there are many future issues to be 
followed after this work as the following: 

1. First of all, it is important to clarify general aspects of jV-fold SUSY in matrix systems, 
as were done in [16, 17] for scalar systems. In the scalar case, there are two significant 
features, namely, the equivalence between A/"-fold SUSY and weak quasi-solvability and the 
equivalence between the conditions (3a) and (3b). In the case of 2 x 2 Hermitian matrix 
systems, however, it does not seem that the conditions (10) and (13) coming from the 
former are equivalent with the conditions (16), (19), (22), and (25) coming from the latter 
although the other conditions are certainly equivalent. That is exactly the reason why we 
considered the both to derive the formula (27). We expect that the general approach [21] 
for the scalar case recently proposed by us would be also efficient for the matrix case. 

2. In the scalar case, the systematic algorithm for constructing an A/"-fold SUSY system [22] 
based on quasi-solvability has shown to be quite effective. Hence, its generalization to 
the matrix case is desirable. It would also enable us to connect directly the possible 
types of matrix A/"-fold SUSY systems with the possible linear spaces of multi-component 
functions preserved by a second-order matrix linear differential operator. For example, it is 
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interesting to see the connection with the quasi-solvable matrix operators constructed from 
the generators of s((2) in Ref. [23]. 

3. Shape invariance is a well-known sufficient condition for solvability [24]. It means in par- 
ticular that it always implies AA-fold SUSY in the scalar case. In fact, some shape-invariant 
scalar potentials were systematically constructed as particular cases of A/"-fold SUSY with 
intermediate Hamiltonians [25, 26]. Recently, several shape-invariant matrix potentials 
were constructed in Refs. [2, 6, 12-14], and we expect that our formulation of jV-fold SUSY 
would be also quite efficient in constructing shape-invariant matrix models. 

4. Extension to more general second-order matrix linear differential operators would 
be possible, e.g., by admitting a non-diagonal second-order operator and by adding a 
matrix-valued first-order operator. In particular, a quantum mechanical matrix model with 
matrix-valued position-dependent mass would be an interesting candidate as a natural 
generalization of Af- fold SUSY in scalar quantum systems with position-dependent mass [27]. 

5. In the scalar case, there are several intimate relations between jV-fold SUSY and A/th- 
order paraSUSY [25, 26, 28, 29]. We expect that we can formulate higher-order paraSUSY in 
quantum mechanical matrix models in a way such that the relations to A/"-fold SUSY in the 
scalar case remain intact in a matrix case. Extension of higher-order A/"-fold paraSUSY [30] 
to matrix systems would be also possible. 
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